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ABSTRACT 

\ 

The  dynamic  stress  intensity  factor  history  for  a  half  plane  crack  in  an 
otherwise  unbounded  elastic  body,  with  the  crack  faces  subjected  to  a 
traction  distribution  consisting  of  a  pair  of  point  loads  that  move  in  a 
direction  perpendicular  to  the  crack  edge,  is  considered.  The  exact 
expression  for  the  mode  I  stress  intensity  factor  as  a  function  of  time  for 
any  point  along  the  crack  edge  is  obtained  by  extending  a  procedure  recently 
introduced  by  Freund  TU  The  method  of  solution  is  based  on  integral 
transforms  methods  and  the  theory  of  analytic  functions  of  a  complex 
variable.  Some  features  of  the  solution  are  discussed  and  graphical  results 
for  various  point  load  speeds  are  presented. 


Availability  Codes 


1.  INTRODUCTION 


A  general  procedure  has  been  introduced  by  Freund  [I]  for  determining 
the  stress  intensity  factor  histories  for  a  class  of  three  dimensional 
elastodynamic  crack  problems.  As  an  illustration  of  the  procedure,  Freund 
studied  a  half  plane  crack  in  an  otherwise  unbounded  elastic  solid,  with  the 
crack  faces  subjected  to  a  pair  of  line  loads  that  are  suddenly  applied 
along  a  line  perpendicular  to  the  crack  front.  Because  the  approach  is 
novel,  its  range  of  applicability  has  not  yet  been  established.  Here,  the 
extension  of  the  procedure  to  situations  with  moving  loads  on  the  crack 
faces  is  considered.  This  distribution  consists  of  a  pair  of  point  loads  that 
suddenly  begins  to  act  at  the  edge  of  the  crack,  and  moves  at  a  constant 
velocity  along  the  crack  faces  in  a  direction  perpendicular  to  the  crack 
edge.  The  corresponding  two  dimensional  problem  was  studied  by  Ang  [2]. 
A  three  dimensional  problem  that  is  related  to  the  one  in  this  paper  is  that 
of  a  point  load  traveling  on  the  surface  of  an  elastic  half  space.  This 
problem  was  analyzed  by  Gakenheimer  and  Miklowitz  [3],  who  considered  all 
point  load  speeds,  i.e.  subsonic,  transonic  and  supersonic.  For  the  purposes 
of  this  paper,  attention  is  restricted  to  subsonic  point  load  speeds,  that  is, 
the  speed  is  less  than  the  characteristic  Rayleigh  wave  speed  of  the 
material.  The  analysis  for  higher  speeds  offers  no  added  mathematical 
difficulty. 

In  section  2,  the  general  formulation  of  the  boundary  value  problem  is 
presented  for  the  three  dimensional  crack  face  tractions  resulting  in  mode 


I  deformation.  Section  3  describes  the  general  approach  to  solving  the 
problem  by  means  of  transfom  methods.  In  section  4,  the  formal  solution 
to  the  particular  traction  distribution  is  obtained  by  means  of  the 
Wiener-Hopf  decomposition  method.  In  section  5  the  dynamic  stress 
intensity  factor  history  is  extracted  by  making  use  of  the  asymptotic 
properties  of  transforms,  the  Cagniard-deHoop  method,  and  the  convolution 
theorem  for  transforms.  A  discussion  of  the  results  is  given  in  section  6. 


2.  GENERAL  FORMULATION 


In  vector  notation  the  Navier  equation  governing  the  displacement  vector 
u  for  an  isotropic  elastic  solid  is  written  as 

0  =  c2jV(V-u)  -  c2sVx(Vxu)  (2.1) 

where  C|  and  cs  are  the  dilational  and  shear  wave  speeds,  respectively.  In 

terms  of  the  Lame  constants  X  and  ji  and  the  mass  density  p,  the  wave 
speeds  are  given  by 

c2,  =  (X+2p)/ p  ,  c2s  =  p/p.  (2.2) 

It  is  also  useful  to  introduce  the  dilatational  and  shear  slownesses  a  and  b, 
where  a=  1/cj  and  b=  l/Cg.  Furthermore,  the  Rayleigh  wave  speed  of  the 

elastic  material  is  denoted  by  cr  and  its  corresponding  slowness  by  r. 

A  standard  approach  when  solving  (2.1)  is  to  introduce  the  displacement 
potentials  <p  and  9  through  the  Helmholtz  decomposition  of  the 
displacement  vector,  i.e. 

u  =  V9  +  Vxy  .  V  9  =  0.  (2.3) 

The  scalar  potential  9  is  called  the  dilatational  potential  and  the  vector 


potential  y  is  the  shear  potential.  The  divergence  free  requirement  on  the 
shear  potential  is  necessary  in  order  to  make  the  decomposition  unique. 
The  advantage  of  this  decomposition  is  that  the  potentials  f  and  ^  satisfy 
the  uncoupled  wave  equations 


<j>  =  c2[V2<p  ,  $  =  c2sV2,f'.  (2.4) 

The  linear  differential  equations  (2.4)  have  the  added  advantage  of  lending 
themselves  to  standard  integral  transform  methods.  The  two  potentials 
are  coupled  through  the  boundary  conditions  that  characterize  the  problem 
to  be  described. 


Figure  1 .  Geometrical  configuration  of  the  elastic  solid. 


Consider  the  elastic  body  containing  a  half  plane  crack  depicted  in 
Figure  1.  A  right  handed  rectangular  coordinate  system  is  introduced  such 
that  the  z  axis  coincides  with  the  crack  front,  and  the  half  plane  crack 
occupies  the  region  y=0,  x<0.  Attention  is  restricted  to  applied  tractions 
of  the  form  Ty  =  ~ cf_(x,z,t)  on  y  =  ±0  where  o_> 0  corresponds  to  a  tensile 
traction.  All  other  components  of  the  imposed  traction  are  zero.  The 
function  cr_  is  prescribed  for  x<0  and  is  extended  so  that  a_  *  0  in  the 

half  range  x  >0.  The  minus  subscript  is  used  to  denote  a  function  that  is 
nonvanishing  in  the  range  x<0.  Likewise,  the  plus  subscript  will  be  used 
to  label  functions  that  are  nonzero  in  the  half  range  x>0,  but  are 
identically  zero  for  x<0.  This  notation  is  useful  in  problems  like  this  one, 
where  the  transforms  of  these  "half  functions"  turn  out  to  be  analytic 
functions  of  the  transform  parameter  in  lower  (minus)  and  upper  (plus)  half 
planes. 

Due  to  the  symmetry  of  the  geometry  and  the  applied  traction,  the 
displacement  fields  satisfy  the  following  conditions: 

ux(x  -y.z.t)  =  ux(x,y,z,t) 

uy(x -y,z,t)  =  -  Uy(x,y,z,t)  (2.5) 

uz(x  ,-y.z.t)  =  uz(x,y,z,t). 

Thus  attention  can  be  restricted  to  the  upper  half  space  y  >  0. 
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Furthermore,  properties  (2.5)  imply  that  aXy(x,0,z.t)  and  cyz(x,0,z,t) 
vanish  for  all  z  and  t,  and  that  Uy(x,0,z,t)  =  u_(x,z,t)  where  u_  represents 
the  unknown  y  component  of  the  displacement  on  the  crack  faces  for  x  <  0 
and  u_  =  0  in  the  half  range  x  >  0.  Hence,  the  complete  set  of  boundary 
conditions  to  be  satisfied  by  the  stress  field  is 

ayyCx.O.z.t)  =  a_(x,z,t)  ♦  a+(x,z,t) 
axu(x.o.z.t)  =  0  (2.6) 

CyZ(X.0,Z.t)  =  0 
Uy(x,0.z,t)  =  u_(x,z,t) 

for  -oo  <  x.z  <  «>  and  0  1  t.  The  function  cf+  represents  the  unknown 
normal  component  of  stress  cfyy  on  x  >  0,  and  <?+  =  0  for  x  <  0. 

The  initial  conditions  are  that  the  material  is  stress  free  and  at  rest 
for  1 1  0.  These  are  expressed  in  terms  of  the  displacement  potentials  by 

<p(x,y,z,0)  =  8t«p(x.y,z,0)  =  y(x,y,z,0)  =  3t9(x,yfz,0)  =  0  (2.7) 


for  y  >  0.  Likewise,  the  boundary  conditions  (2.6)  can  be  replaced  by  their 
corresponding  representations  in  terms  of  9  and  9,  which  are 


xv2<p  ♦  2 M  {  3yy<P  ♦  3y[3z*x  -  3X*Z1  )  =  a.  +  0+ 


(2.8) 


23xy<P  +  3yt3y>fz  -  3zyyl  +  3x[3z^x  ~  =  0 

29yz^  +  ^z^  9z^x  "  9x^z^  +  9y^9x^y  ~  9y^x^  r  0 

3y<P  +  9z^x  "  9x*^z  =  u- 


for  -oo  <  x,z  <  o°  and  0  £  t.  yx,  ^y  and  are  the  rectangular  cartesian 
components  of  4». 


3.  TRANSFORMS 


In  order  to  solve  the  problem  described  by  the  partial  differential 
equations  (2.4),  the  initial  conditions  (2.7),  and  the  boundary  conditions  (2.8), 
use  is  made  of  Laplace  and  Fourier  transforms.  First,  the  one-sided  Laplace 
transform  is  applied  to  the  differential  equations  and  the  boundary 
conditions,  and  the  initial  conditions  are  imposed.  This  Laplace  transform  in 
time  has  the  transform  parameter  s,  and  is  denoted  by  a  superposed  hat.  For 
a  function  <p(x,y,z,t)  this  transform  is  given  by 

r°° 

<P(x.y,z,s)  =  <p(x,y,z,t)  e“st  dt.  (3.1) 

J  o 

For  the  time  being,  s  is  taken  to  be  real  and  large  enough  so  as  to  be  on  the 
right  of  the  abscissa  of  convergence  of  the  integral.  The  dependence  on  z 
and  x  is  then  suppressed  by  applying  a  pair  of  two-sided  Fourier 
transforms.  These  transforms  have  the  transform  parameters  it,s  and  i£s, 
respectively.  They  are  denoted,  respectively,  by  the  upper  case  symbol  for 
the  function  with  a  superposed  hat,  and  the  upper  case  symbol  itself,  to 
wit, 


4>(x,y,C,s) 


-~r  f  <Kx.y.z,s)  eis^z  dz 
n/27T  J  -00 


4>(£,y.C,s)  =  — =■  [  $(x,y,t,s)  e's^x  dx. 

v/ 2rt  J  -oo 


(3.2) 


The  introduction  of  the  s  parameter  into  the  kernel  of  the  transforms  (3.2) 


is  for  algebraic  convenience;  it  takes  care  of  a  change  of  variables  that 
would  otherwise  be  needed  at  a  later  point  in  the  analysis. 

Due  to  the  wave  propagation  character  of  the  solution,  the  strip  of 
analyticity  of  the  Fourier  transforms  can  be  anticipated.  Suppose  that  the 

applied  tractions  are  such  that  a_  vanishes  for  |z|  >  z0.  Then,  at  any 
given  time  t,  the  region  of  causality  is  confined  to  |z|  <  z0  +  tc(.  Thus 

by  considering  the  elementary  wave  field  <p  =  Hit-*- (z+z0)al  + 
H[t-(z-z0)a]  ,  where  H(  )  is  Heaviside’s  unit  function,  and  applying  to  it 
the  t_place  transform,  one  obtains 


0  =  0(e  sOzl-z0)a)  as  |z|  t  oo. 


(3.3) 


This,  in  turn,  implies  that  the  Fourier  transform  <5  converges  for  the  strip 
|lm(C)|<a.  Thus  <J>  defines  an  analytic  function  of  C,  in  the  strip  of 
convergence.  Consequently,  one  can  analytically  continue  <I>  to  values  of  £ 
that  are  not  contained  in  the  strip  of  convergence.  At  this  point  in  the 
analysis  it  is  convenient  to  restrict  C  to  this  strip.  The  identity  property 
of  analytic  functions  [4],  allows  one  to  restrict  C  to  the  portion  of  the 
imaginary  axis  in  the  interval  -a<lm(C)<a,  Re(C)-0.  The  idea  is  to 
perform  the  Wiener-Hopf  factorization  in  the  £  plane  only,  keeping  C 
confined  to  the  strip  of  analyticity  of  6.  At  a  later  point  in  the  analysis  it 
will  be  essential  to  analytically  extend  functions  of  C  away  from  the 
interval  on  the  imaginary  axis  and  outside  the  strip. 


The  domain  of  convergence  of  the  Fourier  transform  in  x  can  also  be 
anticipated.  Suppose  that  the  applied  tractions  are  such  that  along  the 
crack  faces  cL  is  nonzero  for  indefinitely  large  values  of  x  in  the  negative 

direction.  Then,  in  this  region  the  integral  will  converge  provided  that 
lm(£)<0.  On  the  other  hand,  by  definition  the  applied  tractions  do  not 
extend  along  the  positive  x  direction  and  thus,  the  region  of  causality  does 
not  extend  beyond  a  certain  cylindrical  wavefront  ahead  of  the  crack  front. 
To  be  precise,  for  x>0,  the  front  is  centered  at  the  y  axis  and  at  any  given 
time  t,  has  a  radius  z0  +  (x2+ z2),/2.  Thus,  by  considering  the  elementary 
wave  field  9  =  H(t-(z0+  %/x2  +  z2)a)  and  applying  to  it  the  transforms 
(3.1)  and  (3.2)  it  is  found  that  the  final  integral  converges  if 
{lm(£)}2-C2<a2-  Therefore  the  Fourier  integral  $  defines  an  analytic 
function  in  the  strip  -  /C,2  +  a2  <lm(£)<0  in  the  £  plane,  with  C 
restricted  to  -a<lm(0<a,  Re(O  =  0. 

The  class  of  problems  which  is  accessible  by  the  solution  procedure 
outlined  in  the  Introduction  ,  is  one  in  which  a_(x,z,t)  is  restricted  to 
having  a  triple  transform  which  has  the  separable  form. 

—  [  [  a_(x,z,s)e’s(£x+^z)  dzdx  =  —  Z_(£,0  (3.4) 

2 It  J  -oo  J  -oo  sm 

where  m  is  a  real  number  and  2_(£,0  does  not  depend  on  s.  The  reason  for 
this  will  be  more  apparent  in  section  5.  The  requirement  (3.4)  makes  it 
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possible  to  perform  the  final  inversion  of  the  dynamic  stress  intensity 
factor  by  means  of  the  convolution  formula  for  Laplace  transforms. 


The  application  of  the  transforms  (3.1)  and  (3.2)  to  the  partial 
differential  equations  (2.4)  reduces  them  to  ordinary  differential 
equations, 

32<D  a2* 

- s2cx2<l>  =  o  and - s2£2'P  =  o  (3  5) 

3y2  3y2 

where 

<X  =  «(U)  =  -/t?  a2  ,  p  =  g(U)  =  yt.2  ♦  S2  *  b2  (3.6) 

and  a=l/C|  ,  b=l/cs  are  the  dilatational  and  shear  slowneses  introduced 

after  (2.2).  The  complex  £  plane  is  cut  along  s/a2  *  <  |  lm(0  |  <»,  Re(0 

=  0  and  /b2  +  C2  <  |  lm(0  |  <«>,  Re(£)  =  0  so  that  Re(cx)  >  0  and  Re(0)  )>  0 
in  the  cut  plane  for  each  admissible  value  of  £.  With  the  £  plane  cut  in  this 
fashion,  equations  of  the  type  (3.5)  have  the  solutions,  bounded  as  y-»o°  to 
preclude  waves  coming  in  from  remote  regions. 


a(U)  e(U)  qon 

<t>  =  -  e-so<y  ,  S'  =  —  -  e_s^y 


:m+2 


5m+2 


(3.7) 


where  e={B(S.t).C«:.O.D(S.O)  • 


Transforming  the  condition  =0  yields 


SB  -  i£C  +  CD  =  0. 


(3.8) 


It  is  also  necessary  to  transform  the  boundary  conditions  (2.8),  which  will 
yield  four  more  equations.  They  are 


(3.9) 


(b2  +  2C2  +  2C2)A  ♦  2iC£B  -  2i££D  =  *  S+) 

2iCo<A  -  CCB  -  iC3C  +  (a2  +  £2)D  =  0 

2iCo<A  -  (a2  ♦  c2)b  +  icac +  CCD  =  0 

-tx A  -  iCB  +  i£D  =  U_  . 


In  reducing  the  equations  to  this  form  use  is  made  of  the  fact  that  Xcd“2  = 
>i(b2  -  2a2).  Furthermore  we  have  defined  U_  and  as 


5"'  ■  /-oo  roo  . 

U-(S.C)  =  -  u_(x,z,s)e's(^x  +  ^z)  dzdx 

27T  J  -oo  J  -oo 


5"'  r  oo  r  oo  .  /w  >  , 

S+(S.C)  =  -  a_(x,z,s)els^><+^  dzdx 

2 XT  J  -oo  J  -oo 


(3.10) 


The  parameter  s  is  absent,  by  construction,  from  the  five  eqations  (3.8) 
and  (3.9).  There  are  six  unknown  parameters,  four  constants  of  integration 

A.B.C.D  and  two  sectionally  analytic  fuctions  IL  and  Z+,  but  only  five 

equations  to  relate  them.  One  can  solve  for  the  four  constants  of 
integration  in  terms  of  the  two  unknown  sectionally  analytic  functions. 


v 


The  result  is 


A(U )  =  or'(1  -  232b-2)U_(^,0 
B(U)  =  -2iCb-2U_(^0 
c(U)  =  o 

D(U)  =  2i^b-2U_(^,0. 


Substituting  for  these  functions  into  the  first  equation  of  (3.9)  we  obtain 
one  equation  relating  the  two  unknown  functions, 


where 


b2c* 


Rce.ou.a.o  =  s.(U) 


s+(U) 


(3.12) 


R(U)  =  (b2  ♦  2(£2  -  C2)F  "  4(£2  ♦  C2)cx(U)iKUX  (3. 1 3) 

This  is  the  modified  Rayleigh  wave  function,  that  is,  it  corresponds  to  the 
standard  Rayleigh  wave  function,  R(z),  when  C=0.  It  is  a  well  established 
fact  [51  that  the  standard  Rayleigh  wave  function,  in  a  properly  cut  z  plane, 

has  only  two  zeros,  R(iir)=0  where  r  =  l/cr  Thus,  the  modified  Rayleigh 

wave  function  in  the  properly  cut  complex  £  plane  only  has  the  two  zeros 
£  =  ±  i r ,  where  r  =  v/C2+r2. 


Equation  (3.12)  is  a  standard  Wiener-Hopf  equation,  and  the  essence  of 


the  Wiener-Hopf  method  is  to  solve  for  the  two  unknown  sectionally 
analytic  functions  U_  and  £+  from  a  single  equation  (3.12).  The  approach  to 

be  followed  in  section  4  is  actually  due  to  Jones  (6].  Since  £+  is  analytic 

in  the  half  plane  lm(0  >  _  v/T2*"a2~  .  and  U_  is  analytic  in  the  half  plane 

lm(£)<0,  (3.12)  holds  in  the  strip  -i/C2+a2  <£<0.  Furthermore,  C  is 
restricted  to  the  interval  -a<Im(C)<a,  Re(O  =  0,  so  for  a  fixed  value  of  C 
(3.12)  can  be  solved  by  factorization  in  the  £  plane  alone. 


4.  FORMAL  SOLUTION 


At  this  point  it  is  necessary  to  introduce  the  particluiar  applied 
traction  distribution.  The  present  work  is  concerned  with  a  pair  of  point 
loads  that  suddenly  begins  to  act  at  the  edge  of  the  crack  x=y  =  z  =  0  and 
moves  in  the  negative  x  direction,  i.e.  perpendicular  to  the  crack  edge,  at  a 
contant  velocity  v.  Thus,  we  assume  that 


c_(x,z,t)  =  -  P8(z)S(x+vt)  0<v<cr 


(4.1) 


where  S(-)  is  Dirac's  delta  function.  The  amplitude  P  has  physical 
dimensions  of  force  and  P>0  corresponds  to  a  traction  that  tends  to 
separate  the  crack  faces. 


Conical 

headwave 


Dilatational 

wavefront 


Shear 

wavefront 


Figure  2.  Traces  in  the  x-z  plane  of  the  wavefronts  resulting  from  the 
application  of  the  tractions  (4.1). 


Figure  3.  Traces  in  the  y-z  plane  of  the  wavefronts  resulting  from 
the  application  of  the  tractions  (4.1). 


Figure  4.  Traces  in  the  x-y  plane  of  the  wavefronts  resulting  from 
the  application  of  the  tractions  (4.1). 


Figures  2,  3  and  4  indicate  the  leading  wavefronts  that  result  from  the 
application  of  the  tractions  (4.1).  These  wavefronts  have  been  obtained  by 
extrapolation  from  experience  with  two  dimensional  problems  in  hyperbolic 
partial  differential  equations  rather  than  by  solving  for  the  displacement 


potentials.  The  wavefronts  include  spherical  dilatational  and  shear 
wavefronts  centered  at  the  origin  of  coordinates  (Fig.  2).  There  are  also 
two  sets  of  headwaves  that  form  cones  with  vertices  where  the  spherical 
dilatational  front  meets  the  z  axis,  and  extend  to  circles  of  tangency  with 
the  spherical  shear  wavefront.  The  traces  of  these  cones  which  extend  to 
the  crack  faces,  as  well  as  ahead  of  the  crack  edge,  are  shown  in  Figure  2. 
Furthermore,  there  are  conical  headwaves  that  intersect  the  dilatational 
wavefronts  on  the  surfaces  y  =  ±0  for  x<0,  and  extend  to  circles  of 
tangency  with  the  spherical  shear  wavefronts  (Figs.  3  and  4).  The 
heatwaves  arise  since  the  dilatational  waves  alone  cannot  satisfy  the 
traction  free  boundary  conditions  that  exist  at  the  crack  faces. 


Transforming  (4.1)  with  (3.1)  and  (3.2),  it  is  found  that  m=l  and  that 


2-(U) 


iP  1 
2ttv  (C-  ic) 


(4.2) 


where  c=  1/v  is  the  point  load  slowness.  Substituting  for  Z_(£.C)  in  (3.12) 

one  obtains  the  Wiener-Hopf  equation  that  corresponds  to  the  problem  with 
traction  loading  (4.1), 

— -  ♦  I.(i.t).  (4.3) 

b2c<  2ttv  (C-ic) 


Only  some  of  the  steps  involved  in  the  factorization  of  (4.3)  in  the  £ 
plane  will  be  shown  below.  To  begin,  let  ?  =  C2  +  a2  and  note  a  is  some 


constant  in  the  interval  (0,a)  since  C  is  being  held  fixed  in  the  interval 
(-ia.ia).  The  function  o<  can  be  factorized  as, 


<x(U)  =  (t2^V/2  =  [(£M?),/2U(!;-ia),/2]_  .  (4.4) 

The  semicolon  is  used  to  emphasize  the  fact  that  C  is  being  held  fixed  and 
its  influence  in  the  factorization  procedure  is  only  parametric  in  nature. 
The  plus  (minus)  subscript  is  used  to  denote  functions  that  are  analytic  in 
the  upper  (lower)  half  plane  lm(0>  - a"(<0).  The  fact  that  the  two  planes 
overlap  is  no  coincidence.  Indeed,  for  (3.12)  or  (4.3)  to  have  useful 
consequences  there  must  be  a  common  strip  in  which  I+  and  IL  exist.  In 

problems  of  this  sort,  one  hopes  that  cf+  and  u_  are  such  that  the  inversion 

paths  of  their  Fourier  transforms  satisfy  the  above  requirement,  and 
proceeds,  verifying  the  initial  assumption  a  posteriori . 

In  order  to  make  further  progress  one  must  factorize  the  Rayleigh  wave 
function  R(£;£).  For  this  purpose  it  is  most  convenient  to  express  R(£;0  in 
terms  of  a  function  with  neither  zeroes  nor  poles  and  whose  limiting  value 
as  |  £  |  T  oo  is  unity  ([6], pp.  13-15).  Such  a  function  is  easily  constructed 
from  the  asymptotic  properties  of  the  Rayleigh  wave  function  and 
knowledge  of  its  roots.  Specifically,  one  must  make  use  of  the  results 
that 

R($;C)  =  2£2(b2  -  a2)  ♦  0(1)  as  I  $  I  Too 


(4.5) 


and  R(£;0  has  two  symmetric  roots  £  =  ± i  r  ,  where  r  =  /t2  +  r2,  as 
discussed  at  the  end  of  section  3.  The  function  needed  is  easily  seen  to  be 


S(U)  = 


mo 

K (r5  *  £2) 


(4.6) 


where  x=2(b2-a2).  The  function  5  has  the  required  properties  of  being 
nonvanishing  in  the  finite  £  plane  and  being  asymptotic  to  unity  for  large 
values  of  £.  In  the  Appendix  it  is  shown  that  S(£;0  can  be  factorized  as 
S+(£;C)S_(£;C)  where 


S+(£.0  = 


exp1_L  tan  1 
L  tt  J  a 


(b2-2T]2)2 


•qdTi  .  (4.7) 


The  functions  S+(<_,  O  and  S_(£;C)  are  analytic  and  nonzero  in  the  half 
planes  Im(£)>-iF  and  Im(0<5^ 


After  some  straight-forward  manipulations,  (4.3)  can  be  rewritten  as 


fix  U-(£;Q 
b2  G_(£;0 


27tv(£-ic) 


-  E+(£;C)  G+(£;C) 


(4.8) 


where 


G+(U)  = 


(K  ±  ia)1/2 

a±i7)S+(S;0 


(4.9) 


At  this  point  the  factorization  is  almost  complete.  It  remains  to 
remove  the  pole  in  G+(£;C)  at  £  =  ic.  This  is  accomplished  b'j  observing  that 


G+(U) 
K  ~  ic 


G»(frQ  -  G,(ic;Q 
£  -  ic 


G+(ic;Q 
^  -  ic 


(4.10) 


Now  the  factorization  is  complete  and  relationship  (4.8)  may  be 
rewritten  as 


\1K  U_(£;0  iP  G+(iC;0  iP  [G+(£;0- G+(ic;01 


b2  G_(S;C)  2tTv  (^-ic)  2Ttv 


OHc) 


+  G+(£;0£+(£;0.  (4.11) 


The  function  on  the  left  hand  side  of  (4.11)  is  analytic  in  the  half  plane 
lm(O<0.  The  function  on  the  right  hand  side  is  analytic  in  the  half  plane 
lm(0>-iv/t>2  +  a2  •  Since  (4.11)  holds  in  the  strip  where  the  domains  of 
analyticity  of  the  two  functions  overlap,  the  functions  are  each  analytic 
continuations  of  the  other.  Together  they  define  an  entire  function  E(0 
([6],  p.  37)  .  The  function  E(0  is  to  be  determined  from  the  behaviour  of 
the  functions  U_  ,  and  G+  as  |  £  |  T  oo. 
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For  the  class  of  fracture  mechanics  problems  described  in  section  2,  the 
normal  stress  on  the  plane  y=0  is  expected  to  have  the  asymptotic 
behavior 

a+(x,z,t)  ~  kj(t,z)/y/27Tx  asxlO+  (4.12) 

where  kj(t,z)  is  the  dynamic  stress  intensity  factor  history  at  any  point 
along  the  crack  front.  This  history  will  be  determined  explicitly  in  the 
next  section.  In  light  of  (4.12),  <3+(x,z,s)  is  also  expected  to  be  square 

root  singular  as  x  i  0+  for  any  z.  This  result  can  be  used  in  conjuction 
with  the  Abelian  theorem  regarding  asymptotic  behaviour  of  Fourier 
transforms  [61  to  determine  the  asymptotic  properties  of  £+(£;C),  namely, 

lim  /2£/s  e~ iTt/4  £+(£;0  =  lim  /x"a+(x,z,s)  .  (4.13) 

x  1  o+ 

Since  the  right  hand  side  of  (4.13)  is  equal,  by  (4.12),  to  a  function 
parametric  in  s  and  C.  £♦(£;(,)  =  0(£~,/2)  as  £too.  Furthermore,  u_  is 

expected  to  vanish  as  xTO"  and  |  G±(^;C) |  =  0(|£|-,/2)  as  |£|Too. 

Therefore,  both  sides  of  (4.11)  vanish  as  |£|  Too.  That  is,  E(0  -»  0  as 
|^|  Too.  Fr  im  the  preceding  discussion  it  follows  that  E(£)  is  a  bounded 
function  in  the  entire  £  plane.  According  to  Liouville's  theorem  a  bounded 
entire  function  is  a  constant  [41,  and  therefore  E(£)  must  be  identically 


a 


zero,  i.e.  E(O-0.  Equation  (4. 11)  can  now  be  solved  Tor  I+  and  U_,  which 


are  found  to  be 


„  x  P  t>2  G+(ic;0G_(^;0 

U_(U)  = - - - 

2Ttiv  pK  £  -  ic 


£♦(£*)  =  - 


1  rG+(ic;0 


(4.14) 


27tiv  (£-ic)  LG+(£;C) 


This  completes  the  formal  solution  of  the  problem.  In  theory,  the  exact 
solution  to  the  problem  is  found  by  taking  the  triple  inverse  transforms  of 
(4.14). 
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5.  INVERSION 


Unfortunately,  an  exact  inversion  of  Z+  and  U_  as  given  by  (4.14)  is  not 
evident.  This  section  is  concerned  with  the  determination  of  the  stress 
intensity  factor  history,  which  is  given  by  the  asymptotic  inversion  of  £+ 
with  respect  to  £.  As  mentioned  in  section  4,  the  asymtotic  behaviour  of 
I+(£;C)  as  |^|  Too  is  intimately  related  to  the  behaviour  of  the  double 

transform  Z+(x,0.  of  the  normal  stress  a+(x,z,t)  as  x  i0+.  Jt  follows  from 
(4.12)  that  the  Laplace  transform  of  the  dynamic  stress  intensity  factor 
history  k jCt.z)  is  simply 


kj(s.z)  =  lim  v/2xrx  a+(x,z,s).  (5.1) 

xiO+ 

Thus,  the  double  transform  of  k|(t,z)  is  found  from  (4.13)  and  (4.14)  to  be 

P  v/c+(t,2  +  a2)1/2 

Kj(s.C)  r - v/tt/s  - .  (5.2) 

Ttv  [c+(C2  +  r2)1/2lS+(ic,0 

By  construction,  the  function  Kj(s.C)  is  analytic  on  the  strip  -ia<C<ia. 
Since  the  Wiener-Hopf  factorization  has  been  completed  and  the  limiting 
process  that  resulted  in  K|(s,C)  has  suppressed  all  dependence  on  the 

variable  £,  there  is  no  longer  any  need  to  restrict  C  to  this  strip.  For  the 
purposes  of  inverting  this  function  it  will  be  advantageous  to  extend  the 
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definition  of  the  function  to  the  entire  C  plane.  It  is  this  extended 
function  which  is  referred  to  as  the  Fourier  transform  of  k|(s,z).  This 

function  can  be  made  single  valued  by  cutting  the  C  plane  along 
a<  ]lm(0]  <<*>,  Re(O  =  0  and  along  r<  |lm(0|  <0°.  Re(t)  =  0.  This  ensures 
that  the  square  roots  have  positive  real  parts.  The  function  S+(ic,C)  is 

analytic  over  the  entire  C  plane  by  virtue  of  a  theorem  pertaining  to 
analytic  functions  defined  by  integrals  ([4},  p.  92;  [6],  pp.  I  1-12). 

The  inverse  Fourier  transform  of  (5.2)  is 

S  r  . 

k|(s.z)  =  —r  K|(s.C)e~ls('zdt  (5.3) 

where  Co  is  3  real  number  between  -a  and  a.  it  is  hoped  that  the  final 
Laplace  inversion  of  k](s.z)  can  be  done  by  means  of  the  convolution 

formula.  Therefore,  it  will  be  advantageous  to  cast  (5.3)  in  the  form  of  a 
one  sided  Laplace  transform  so  that  the  inversion  can  be  performed  by 
inspection.  The  required  transformation  can  be  achieved  by  suitably 
deforming  the  Fourier  inversion  path  into  a  branch  line  integral;  this  is  a 
trivial  case  of  the  Cagniard-deHoop  technique.  Consider  the  case  of  z>0. 
For  (5.3)  to  be  a  convergent  integral  when  z>0,  it  is  necessary  that  the 
inversion  contour  lie  in  the  lower  half  plane  lm(C)<0.  Since  the  finite  C 
plane  does  not  contain  any  singularities  aside  from  the  branch  points  at  C 
=  Jiia  and  C  =  ± ir ,  the  original  inversion  path  can  be  deformed  into  a  new 
contour  as  depicted  in  Figure  5.  This  new  contour  consists  of  two  quarter 


A 

Figure  5.  The  complex  %  plane  showing  the  singularities  of  K|(s,$)  and 
the  integration  path  for  the  evaluation  of  (5.3). 


circles  r,  and  r2  in  the  lower  naif  plane,  and  a  branch  line  path  running 
upward  along  the  left  hand  side  of  the  branch  cut  from  -i<»  up  to  the 
branch  point  at  C  =  -ia,  around  the  branch  point,  and  finally  running  down 
the  right  hand  side  of  the  branch  towards  That  the  original  integral 
in  (5.3)  is  equivalent  to  the  integral  of  the  same  argument  along  the  new 
path  follows  from  Cauchy's  integral  theorem.  This  theorem  states  that  the 
contour  integral  of  a  function  which  is  analytic  on  and  inside  the  contour 
is  identically  zero  HI.  The  contribution  of  the  integral  along  the  arcs  at 
infinity  Tj  and  r2  as  |C|  to°  vanishes  by  Jordan's  Lemma  HI.  Thus  (5.3) 
reduces  to  a  branch  line  integral.  Exploiting  the  fact  that  Kj(s,-C)  =  K(s,0. 

where  the  bar  denotes  complex  conjugation,  the  branch  line  integral  can  be 
expressed,  after  an  elementary  change  of  variables,  as 


J  e_STl  !m{F(+0+iTi/z)}dT| 


[  c  +  ( C2  +  r2 )  ,/2]S+(  ic,0 


In  (5.4),  F(C)  is  evaluated  on  the  right  hand  side  of  the  branch  cut.  The 
s-multiplied  Laplace  transform  in  (5.4)  can  be  expressed  as  the  Laplace 
transform  of  a  derivative,  i.e. 


kj(s,z) 


/2  P 


oo 

z-i  s~\/2  f  e-ST|  — lm{F(+0+i'f|/z)}d'q  (5.5) 

J _  Rti 


There  is  no  endpoint  contribution  from  T)=az  in  (5.5),  because  F(+0+ia)  is 
a  real  quantity.  The  inversion  of  (5.5)  is  now  obvious  because  kj(s,z)  is 

seen  to  be  the  product  of  two  transforms.  Therefore  kj(t,z)  is  the 
convolution  integral  of  the  inverses  of  the  two  transforms, 

/2P  rt/z  q 

k](t.z)  =  — —  z’1  —  lm{F(+0+iO)  .dC  H(t-az)  (5.6) 

tt3/2v  Ja  sc  y  t  -  Cz 


for  z>0,  and  k|(t,-z)  =  kj(t.z).  The  expression  for  kj  in  (5.6)  apparently 

cannot  be  reduced  further  in  terms  of  elementary  functions.  Some  of  the 
properties  of  the  real  integral  (5.6)  along  with  an  interpretation  of  the 
results  are  discussed  in  the  next  section. 


6.  CONCLUSION 


Even  though  the  dynamic  stress  intensity  factor  history  (5.6)  cannot  be 
evaluated  in  terms  of  elementary  functions  it  can  be  evaluated  numerically 
and  some  of  its  salient  features  can  be  obtained  analytically.  Figures  6,  7, 
8  and  9  show  the  results  of  the  numerical  integration  of  the  integral  (5.6) 
for  the  values  of  the  ratio  of  the  point  load  velocity  to  the  Rayleigh  wave 

speed  of  0.8,  0.6,  0.4,  and  0.2  for  a  Poisson  ratio  of  v  =  0.3  (cj/cr  =  2.02). 

The  time  scale  has  been  non-dimensional ized  so  that  r-  1  corresponds  to 
the  arrival  of  the  dilatational  wave  at  the  observation  position  z  along  the 
crack  front.  The  dynamic  stress  intensity  factor  has  been  normalized  by 
premultiplyng  (5.6)  by  P_1(Tfz)3/2. 

Following  the  sudden  application  of  the  point  loads,  a  point  z  along  the 
crack  edge  is  at  rest  until  the  arrival  of  the  dilatational  front.  This  front 
is  compressive  in  nature  and  the  crack  faces  respond  to  it  by  initially 
trying  to  close  together.  This  is  reflected  by  the  stress  intensity  factor 
being  negative  initially.  The  initial  jump  in  the  dynamic  stress  intensity 
factor,  which  can  be  seen  in  the  graphs,  is  a  verifiable  feature  of  the 
solution  of  this  particular  three  dimensional  loading  distribution.  For  the 
case  of  a  pair  of  line  loads  acting  perpendicular  to  the  crack  edge,  i.e  the 
convolution  of  the  traction  distribution  (4.1)  for  v  =  0,  Freund  Ill  found  that 
the  dynamic  stress  intensity  factor  started  from  zero  and  gradually  became 
negative.  By  taking  the  limit  as  t/z  i  a+  in  (5.6)  the  jump  is  found  to  be 
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lim  ki(t,z)(7rz)3/2  P~'  = 
t/z  l  a+ 


A  [c  + 


|S+(ic,0++ia) 


(6.1) 


The  dynamic  stress  intensity  factor  becomes  increasingly  negative 
until  the  arrival  of  the  Rayleigh  waves  at  z-2.02  when  it  becomes 
logarithmically  singular.  Between  r  =  1  and  r  =  2.02,  kj(t,z)  also  exhibits  a 
disturbance  coinciding  with  the  arrival  of  the  shear  front  at  r=  1.88.  This 
mild  discontinuity  in  the  slope  of  kj(t,z)  is  due  to  a  change  in  the  form  of 

the  function  S+  as  the  branch  point  located  at  the  shear  wave  slowness  is 

crossed.  After  the  passage  of  the  Rayleigh  waves  the  crack  faces  begin  to 
open,  the  stress  intensity  factor  increases  until  reaching  a  maximum,  and 

thereafter  decaying  very  gradually  towards  its  limiting  value  kj(°°,z)  =  0. 


This  completes  the  analysis  of  the  three  dimensional  stress  intensity 
factor  history  for  the  case  of  a  pair  of  moving  point  toads  on  the  faces  of 
a  crack.  The  solution  (5.6)  to  the  problem  described  by  the  traction 
distribution  (4.1)  is  the  fundamental  solution  to  the  class  of  problems 
involving  traction  distributions  moving  perpendicular  to  the  crack  edge.  In 
this  paper  the  range  of  applicability  of  the  procedure  introduced  by  Freund 
[I]  has  been  successfully  extended  to  this  class  of  problems.  Other 
situations,  such  as  moving  loads  along  a  direction  inclined  to  the  crack 
edge,  can  be  examined  following  the  same  methodology  used  in  this  paper. 
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re  8.  The  normalized  stress  intensity  factor  kj(t.z)  (ttzP^  p~1 
versus  t/az  for  the  case  of  c/r  =  0.4. 


versus  t/az  for  the  cose 
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of  c/r  =  0.2. 


APPENDIX 


In  this  appendix  the  factorization  of  the  function  S(£;0  as  defined  by 
(4.6)  is  carried  out.  As  mentioned  in  section  4,  the  function  S(£;C)  has  no 
zeros  or  poles  in  the  £  plane  cut  along  F<|lm(0|<b  ,  Re(O  =  0. 
Furthermore,  the  function  S(£;0  has  the  desired  property  of  being 
asymptotic  to  unity  as  |  £  |  T  oo.  These  properties  enable  one  to  factorize 

S(£;C)  into  the  product  S+(£;t)S_(£;0  by  a  well  known  decomposition 

technique  ([61,  pp.  15-17).  Specifically,  one  considers  the  Cauchy  integral 
representation  for  lnS(£;C),  that  is,  one  represents  In  S(£;0  by  the 
following  contour  integral 

lnS(£;0=  — ■ —  <j>  ln5(tf;t)_£l2S__  .  ( A 1 ) 

2ni  J r  V-Z 

The  counterclockwise  contour  r  is  depicted  in  Figure  10.  By  expanding  the 
contour  r  to  infinity,  making  use  of  the  asymptotic  behaviour  of  S($;0.  an 
equivalent  representation  of  (Al)  is  obtained.  The  integral  (Al)  around  the 
contour  r  at  infinity  is  equivalent  to  two  clockwise  contours  r+  and  r_ 
around  the  branch  cuts  as  shown  in  Figure  10,  so  that 

lnS(i;t)  =  In  S,(£;0  *  In  S-(M)  =  — —  I  lnS(y;t)  <U_  .  (A2) 

2iri  J  r.  *  r.  y-t 


Thus,  the  functions  S+(£;C)  and  S_(£;C)  are  given  by 


where  2f=lm(0-  Along  the  left  hand  side  of  r+,  R(£;0  is  given  by  the 
complex  conjugate  of  (A4).  From  (4.6) 


lnS(£;0  =  lnR(£;0  -  ln[ 2 (b» - a^X r  +  t2)) 


(A5) 


and  therefore  S+(£;C)  can  be  expressed  as, 


In  S+(£;C)  =  -1-  f  lnR(+0+  iff)  dtf  -  —  f  1  lnR(-0+iy)  d 7S  (A6) 

2rt  J-a  ftf-?  2jt  J-a  fa -7 


=  —  [  ^  In  [r(+0+  i2f)/R(-0+  i#)1  dff 
2tt  j  -  a  L  J  Tx^i 


(A7) 


i  r  b 


7t  J  a 


tan-1 


(b^H2-22f2)2 


d2f  .  (A8) 


The  function  S_(£;0  is  treated  similarly.  One  further  change  of  variables 
T[  =  (tf2-C2)1/2  reduces  (A8)  to  the  form  (4.7). 
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